Stationary states of Jaynes-Cummings model with atòmic center-of-mass quantum 
motion: direct comparison of standing-wave and counterpropagating-waves cases. 
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The eigenstate problem of the Jaynes-Cummings model on the basis of complete Hamiltonian, 
including the center-of -mass kinetic energy operator, is treated. The energy spectrum and wave 
functions in standing-wave (SW)- and counterpropagating waves (CPW)- cases are calculated and 
compared with each other. It is shown that in CPW-case i) the atòmic momentum distribution is 
, asymmetric and somewhat narrower in general; ii) the concept of quasimomentum is not applicable 

and instead the ordinary momentum concerns the problem; iii) atòmic and photonic state distribu- 
tions are self-consistent, and, in consequence iiii) mean number of photons in the counterpropagating 
O ■ traveling waves and mean atòmic momentum match. Explícit analytic expressions for energy eigen- 

' vàlues and eigenfunctions are found in Tavis -Cummings-type approximation [Phys. Rev. 170, 

U—\ , 379(1968)] and is pointed, that it implies only the bounded-like states for atòmic center-of-mass 

\q • motion. It is also shown that if the recoil energy is taken into account, the Doppleron resonance is 

| split into two branches, one of which diverges to Bragg-like resonance in the high-order range. 
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> ! I. INTRODUCTION 

The key scheme of cavity quantum electrodynamics (QED) [Q, modern atom òptics and interferometry Q is the 
' resonant interaction of an atom with cavity fields, created as standing or counterpropagating waves. The first type 
t-H , (with respective quantization on the SW basis) is attained in two plane-parallel-mirror cavities (in short, cavity), 
while the second type (with respective CPW-quantization) is attained in three or more mirror cavities (ring cavities). 
It is well known, in addition, that for classical picture of fields, these two representations are equivalent in the sense, 
■ that the SW can always be presented as a superposition of two counterpropagating travelling waves. In the quantum 
O theory, nevertheless, they are divers [||, including the Hamiltonians of interaction. Hence it needs to be ascertained: 
, 1 1 . are the diversities only quantitative or qualitative too, and how much are they for these or that circumstances and 
processes. Note, that it presents not only an acadèmic interest, since in the microcavities (V ^ 10 _3 cm 3 ) several 
dozens of photons can induce strong optical nonlinearities for dipole-allowed optical transitions. 

As far as we know, elucidation of physical aspects of mentioned incongruity had been done in paper 0] , taking as an 
example the process of near-resonant coherent diffraction of atòmic matter-waves by a space-periodic làser radiation 
field, known as near-resonant Kapitza-Dirac diffraction Q . By means of numerical solution of the master equations 
for atòmic probability amplitudes, the existence of evident differences between SW and CPW diffraction patterns was 
shown , in strictly quantum domain of cavity fields (the mean number of photons in cavities was chosen one). The 
physical reason of the difference also was presented. It is the behavior of momentum conservation law in relevant 
"atom + field " systems. This general principle ceases to apply in SW-case, because the atom changes its own 
momentum due to the interaction, while the SW-field does not. In contrary, in CPW-case both, atòmic and field 
subsystems, change the momenta equally and in opposite directions, thereby conserving the systems total momentum. 
Already on the base of this sole difference can be anticipated a " soffcr" diffraction for CPW-case with a respectively 
narrower distribution of states in the momentum space. 

In this paper are considered the stationary states of the system "two-level atom + quantized field" in both SW and 
CPW-cases and the qualitative and quantitative differences between them are carried out in the range of quantum 
òptics (the number of photons should be less than several tens). To this end the simplest theoretical model, the 
Jaynes-Cummings model (JCM) including the quantized atòmic center-of-mass motion [Q, will be used. The 
spontaneous emission and other incoherent processes, as usual for this model, are omitted. 

To make possible the direct comparison of these two concepts we, at first, represent detailed theory for stationary 
states of the systems at hand. Our representation for SW-case somewhat coincides with the picture in ||. The 
preliminary analysis of the states and the comparison of the SW and CPW-cases are based on numerical solutions 
of the exact, as well as of the off- resonant approximate equations of atòmic probability amplitudes in momentum 
space. An additional approximate analytical solution of the problem allows us to find out the parameters (terms), 
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determining the size of nonequivalency between the SW and CPW-spectrums and eigenfunctions in explícit form. 
The latter solutions pertain, however,only to bounded (wells' inner) atòmic states. 

The paper is organized as follows. Section II represents the bàsic model and the exact set of equations. Here 
we also give some numerical simulations of the problem, implemcntcd for the sodium atom with 3S]/a — 3P 3 / 2 main 
transition. Section III represents the off-rcsonance and mentioned analytical approximations. In Section IV we 
examine the Bragg and Doppleron resonances, including the recoil energy terms. In Section V we consider the 
mechanisms responsible for the formation of stationary states from free definite-momentum ones. Conclusions and 
some short remarks about the subject are given in Section VI. 



II. BASIC MODEL THEORY 



We consider a two level atom of mass M and optical transition frequency lüq, interacting with a quantized plane 
monochromatic field of frequency u>. We will restrict ourselves with two important cases of standing and coun- 
terpropagating waves. The longitudinal components p x ,p y of the atòmic center-of-mass (cm.) momentum remain 
unchanged for perfectly plane fields, so that only the transverse atòmic momentum p z = p needs to be considered. 
Since the incoherent processes are not included, the system is governed by the Schrodinger equation. 

Let us first consider the CPW-case. The Hamiltonian in dipole and rotating-wave approximations is given by Q 

H = TTjrjP 2 + \ í 1 + °3) fi^o + ^ (afai + cà,a 2 ) 



+ 3 ^ 2 hílo [er + (ai exp(ikz) + a 2 exp(— ikz)) + er_ (a± exp(ifcz) + exp(— i/cz))] , (1) 

where P = —ihd/dz is atòmic momentum operator (1D), 03, a+ and o— are usual pseudospin atòmic (Pauli) operators, 
a-i and (i = 1, 2) are annihilation and creation operators for two running wave modes, k — ui/c. The factor ÍIq in 
coupling constant is the vacuum Rabi frequency, and is connected with the atòmic transition dipole moment d by 
relation ü,q/2 3 ^ 2 — —d(2TTu>/hV), where V is the normalization (cavity) volume. The coherent binding between the 
internal and cm. variables is by the Rabi frequency f2 , too. 

The system under consideration has four degrees of freedom: two per atòmic internal and 1D cm. motions, and 
two per counterpropagating travelling waves. Three operators, forming together with the Hamiltonian the complete 
set of mutually commutative operators, are known and are the following: the excitation number operator 

N =^{l + a 3 ) + a+ ai +a+a 2 , (2) 

the total momentum operator 

P = p+ hk (a^ai — a^a 2 ) , (3) 

and the operator 

which combines translational p and dipole inversing 03 operators . Therefore H, N, P and T have compatible 
eigenvalues and a common system of eigenfunctions. Only these states, being the bàsic, should be considered 
thereafter in this paper. Any other state, of course in principle, can be presented via these eigenstates. Note, that 
sometimes the operator / = exp(i27rp/?ïfc) is used in related problems. However, / = T 2 and can not be added into 
the set of mentioned operators or replace the operator T there. 

Out of a desire to better understand the future solutions and their connection with free-state vàlues, which present 
a definite interest too, it may be worth to emphasize here that all three N, P and T operators don't contain the 
interaction parameter ÇIq. They are interaction independent. Their eigenvalues, denoted by TV, P and T respectively, 
also are interaction independent. Therefore, any intermediate or final state, created from the initial one due to 
interaction should possess the same vàlues of N, P and T as the free initial state. These ones, of course, don't need 
to be stationary in general, but may be such, if the interaction is adiabatic || . A sufficient condition for adiabatic 
following is the exceeding of the interaction switching on and switching off times (t sw ), the inverse quantity of the 
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system's smallest characteristic frequency (íí^in)- ^ n the system at hand, as a such frequency appears the one-photon 
recoil frequency íl r = hk 2 /2M. So for the interaction times essentially exceeding the recoil time íl" 1 , the atòmic state 
evaluation is adiabatic with respect to both, internal and cm., motions of the atom |ïc|] . Failure of the adiabatic 
condition implies, respectively, that the population leaks into the immediate neighboring states, first of all into the 
translational states. Applying the general remarks to the system under consideration, we arrive to the possibility 
to treat the stationary states as attained adiabatically from some free initial state and then replace the set of N, P 
and T quantities by a new set of three related parameters of this initial state: atòmic moment um p and the photon 
numbers niand n% in counterpropagating waves. One must keep in mind, however, that the mentioned replacement, 
which will be embodied soon, is substantial for adiabatically treated stationary states and is not obligatory even in 
this case. 

Dirac notations \l) 1 and \m) 2 are used for quantized travelling-wave states (number states). Atòmic internal ground 

and excited states are specified by ( ^ J and í ^ J matrixes respectively, the cm. states- by the function exp(ipz / K) . 

Hilbert space of the " atom+field" system is the direct sum of Hilbert spaces of atòmic cm. and internal motions, as 
well as of the field modes. Thereby, the nondegenerate wave function of the free system can be written as 

*/ ree = Aipjexp I — j \l) 1 \m) 2 , (5) 

where A is the normalization constant, j = g, e; <p gte — í ^ ^ , í J ^ ; — oo < p < oo and i, m = 0, 1, 2, • • • in general. 

The seeking wave functions of the interacting system can be written as a superposition of (|^)-type terms. 
The general form of a wave function with a definite N (N*S? — N^) is 

N 



*=(?)£/ UrM^^N-r^i^dq 



N-l 



) E Mí.ílHilW-l-^exppf dç, (6) 



r=0 

where a r and b r are arbitrary coefficients yct. Admissible solutions also should be normalized according to \ "f) 
= 1. Note, that the excitation number N can be represented by means of the numbers n\ and (see the remark 
above ): 

„ _ ( «i + n 2 for free ground level atom, \ . , 

1 ni + n-2 + 1 for free excited level atom. J ' ^ ' 

validity of which can be directly checked with the help of wave function (|^) . 

For the next step we demand from ^ to be the eigenfunction of total momentum operator P: 

= P*. (8) 

Inserting (|) and (§) into (|) we arrive to equations 

a r (q, t) = a r (t)S {q + (2r - N)hk - P) , (9a) 

b r (q, t) = b r (t)S (q + (2r + 1 - N)hk - P) , (9b) 

where S(x) is the Dirac delta - function and a r (t) and b r (t) are the new coefficients. The eigenvalue P can be written 
as : 

P = p+{n 1 ^n 2 )hk 1 (10) 

where p is chosen the same for both, ground and excited atòmic levels. Integration over q now yields to the following 
form: 
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N 



*=( 1 )T, a r( t )\ r )i\ N - r )2^(kP-(N~2r)hk) : . 



(11) 



JV-1 



1 ) Y J K{t)\r) l \N-l-r) 2 c^( l -{P-{N-l-2r)hk) 



This form is already eigen-one for the operator T 



T = -exp(-(P-^)l=-exp^-p 



) I = - exp ( 



(12) 
(13) 



Hence, the expression ( |ïl·| ) is the eigenfunction of operators AT, P and T . Using(Q), ([ÏÏ]) and the second equation 
in (13), the latter three can be replaced by the initial vàlues n\,ri2 and p. 

As it was previously stated, we are concentrating on the stationary states, where 



i iEt\ , , , , / iEt\ 
a r (t) = a r exp [ — I , b r (t) = b r exp I — I , 



(14) 



E being energy of the system. 

After Standard transformations the Schrodinger equation leads to following set of tridiagonal recurrence algebraic 
equations: 



E 



Nfvjj - (P - (N - 2r) hkf /2M a r = (hfl /2 3/2 ) Vrb r -i + \/N - rb r 



(15a) 



E + hs- Nhoj - (P - (N - 1 - 2r) hkf /2M b r = (hn /2 3/2 ^ V7Tïa r+ i + y/N 



(15b) 



where r = 0, 1, 2, ■ ■ -, N and e = lo — lüq is the atom-field frequency detuning. 

Equations ( |15a| ), ( |15b| ), as well as equations (18a) and ( |18b| ) for the SW-case, are the basis of the remainder of 
this paper. They determine the coefScients (probability amplitudes ) a r , b r for ground and excited internal energy 
levels and the per mitte d vàlue s of energy E. 

The set of Eqs. (|Ï5a|) , ( pb| ) has 2N + 1 solutions E„(P, N), {a r {v, P, N), b r {v, P,N)\u = 0, 1, 2, 2A^}. (N + 1) 
of them may be interpreted as created adiabatically from initial states with a ground-level atom, and the other N ,as 
created respectively with an excited-state atom. For simplicity we suppress N, P and v in probability amplitude and 
energy notations. 

Rather than pursuing the governing equations in this CPW-case we now turn to another case, to SW-case. The 
respective Hamiltonian is H] 



H = T^J^P 2 + 2 ^ ^ UjJ ° ^ lüja+a + -^Klo{<J+a + <j-a + ){exp(ikz) + exp(— ikz)], 



(16) 



where a and a + are now SW-mode annihilation and creation operators, obeying the same-type commutation relation 
[a, a + ] = 1 as the running modes. Note, that the coupling constant in this case is smaller by factor 2 1 / 2 than in the 
former case. 

The system "atom + SW" has three degrees of freedom, by one less than the system with CPW. The reason is, 
that the mo mcnt um P = p,and is not commutative with the Hamiltonian (|Ï6|). Proceeding the same way as for ( |lT 
and (15a), ( 15b ), we arrive to 



*= ( 1 ) M E a l (t)exp( ^(p-2lhk) 



l= — oo 



(17) 



l) 1^-1) E b l {t)ex V ( l -(p-{2l + \)hk) 
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for wave function, and 



E - Nhiu - (p - 2lhk) z /2M oí = (fiílo/4) V N [6;-i + h 



(18a) 



E + he- Nhaj - (p - (2/ + 1) /2M 6/ = (ftfi /4) ViV [aj+i + aj] (18b) 

for probability amplitudes and energies. Here l — 0, ±1, ±2, ... , and the other notations are obvious. For adiabatically 
attained states 

j. _ ( n for free ground level atom, 1 , . 

I ri+1 for free excited level atom. [ ' ^ ' 



Besides the apparent similarities, the bàsic sets of equations, (15a), (15b) and ( |l8a| ), ( |l8b| ), also have bàsic dif- 



ferences. While the first set is finite-size, the second one is infinite (even for N = 1, one-photon case). The next 
difference is in right-hand side coefhcients. They are varied (r-dependent) for the first, but are constant for the 
second set. And the final difference, which already is conceptual and deserves a peculiar attention, is that the total 
momentum P in CPW-case, in contrast to its counterpart p in SW case, cannot be treated as a quasimomentum. 
This as serti on immediately follows from the fact, that the replacing P — > P + s2hk 1 with the relabeling r — * r — s in 



(15a), ( 15b ), s being an integer, conserves the left-hand side coefficients invariant, while the right-hand side coefhcients 
get new vàlues, leading thereby to other vàlues for eigenenergies and probability amplitudes. Reminding that the 
physical reason of this difference is the participation of the momentum conservation law, we arrive at a conclusion 
that the admission of the quasimomentum concept in space — periòdic systems is conditioned by the failure of 
momentum conservation. There is no concept of quasimomentum in really closed systems conserving the total 
momentum. 

The main goal of the remainder of this section will be to illustrate the behavior of both systems via the results of 
numerical calculations. While the solutions of the CPW equations have been got exactly, to SW equations has been 
applicd suitably truncated matrix diagonalization mcthod. Solutions depend, except the number N of excitations 
(photons), upon the relative size of four parameters: the resonance detuning e, vacuum Rabi-frequency Q, transverse 
kinetic energy detuning Pk/2M and photon-induced kinetic detuning, which is the recoil energy E r = (2hk) 2 /2M 
in h units . It is convenient to introduce dimensionless parameters, scaling the energies by the recoil energy E r and 
all momenta - by one-photon-reemission recoil momentum 2hk. Then, as dimensionless parameter of interaction for 
both sets of equations appears 

£ = hn /2 3/2 E r . (20) 

The strength of coupling is essential already at £ ~ 1. Let us for illustration take the sodium's atom with D2 
(3Si/2 — 3P3/2) main transition and estimate the volume V of microcavity for which this is so, more definitely: 
^Na — 10. Using the explicit expressions for Qq and E r we arrive at 



( , 3/ |e| MX 2 

t = V o~i7 9 2 fe ' 20a ) 
V 2mV 2n^h 

where / and A are the oscillator strength and the wavelength of optical transition, m is the mass of electron. Inserting 
the vàlues (A = 589 * lCT 7 cm, M = 3.8 * 10~ 23 g, / = 0.655) we get V = 10~ 3 cm 3 . Note, that for cesium atom, with 
the analogous transition 6Si/ 2 — 6P3/2 and for the same volume of microcavity, the coupling strength is more than 
one order stronger, Çcs ~ 135. 

Figures la and lb show the energy spectrum as a function of momentum for CPW and SW-cases respectively. In 
both cases N — 6, V — 10~ 3 cm 3 , he/E r = 250, and the atom of Na is chosen. The comparison of figures clearly 
indicates the mentioned diffcrcnccs for energy spectrums: first, the number of energy branches is finite (2A^ + 1 = 13 
in chosen case) for CPW-case, whilst it is infinite for SW-case (are presented only some lower branches); second, the 
CPW-spectrum is not repeated as 2hk -wide zones, pointing out the failure of the quasimomentum concept in this 
case. The Figure || , plotted for the case N — 12, shows, as it was expected, approaching of the CPW-picture to 
the SW one. The flattening of lower curves signifies the transition of these states into the bounded-like states, for 
increased number of external photons. One can see, that the gaps between the energy zones are essentially wider in 
CPW-case. Also is seen the well known regularity, that the widths of the gaps between the n th and (n + l)*' 1 bands 
diminish with number n |ïï| . 

Probability amplitude distributions are shown in Figures 3a, 3b and 4a, 4b, respectively for N — 6 and A^ = 12 
cases. As is seen from comparison of the graphs, in the first pair ( A^ = 6 full quantum range ) there is an apparent 
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discrepancy between CPW and SW-cases. In the second pair of graphs (N — 12) it is suppressed, but yet stays 
visible, even essential. The latter graphs also display the difference of momentum distributions for bounded-like 
(lower-laying levels) and free-like (upper laying levels) states. 

The momentum distribution in bounded-like stationary states has a two-peak form. The peaks retire from each 
other for increasing energies, gradually being broader and losing in height. In free-like states , in contrast, the distri- 
butions are centered at the nonperturbed vàlues (i.e. at l — in SW-case). To avoid the eventual misunderstanding 
of depicted graphs, it needs to be noted that the momentum axes (abscissa) in these figures, as well as in the following 
ones for probability amplitude distributions, really contains only discrete vàlues. Only for convenience of exhibition 
the neighboring points in graphs have been merged via the straight lines 

III. OFF-RESONANCE APPROXIMATION AND TAVIS-CUMMINGS-TYPE ANALYTIC SOLUTION 



In this section we examine the case, where the resonance detuning e is much larger than any other above men- 
tioned frequeney; íïo 1 P(p)k/2M and E r /h. This case u suall y is r egard ed as adiabatic elimination, or off-resonance 
approximation p2|. In zero order for e we obtain from (15b) and (18b) the equations 



(n /2 3/2 e) V7Tïa r+ i + \/ÏV-~ 



(21) 



r = 0,1,..., N 



k ~ (VNÍl /4e) [a l+1 + 0j ] 



(22) 



l = 0,±1,±2, ... respectively. This approximation, widely used in problems without the kinetic energy operator in 
Hamiltonian, calls for some comments here. The matter is, that for instance in SW-case, the term (p— (2l + l)hk) 2 /2M 
already can exceed the any preliminarily given value he. But it has to be taken into account, that the probability 
amplitudes, rather out of two-peaked range of distribution, are extremely small and can be put out of consideration. 
Then we arrive at the condition 



£ » 



(23) 



for implemcntation of the approximation, the same as the one without the kinetic energy operator: detuning of the 
resonance must be much greater than the o ptica lly in duced inhomogeneous width of the transition. 

By means of ( ^l·| ) and ( p2[ ) the Eqs. (15a) and (18a) are being transformed into the recurrence equations only among 
the ground-levcl atòmic amplitudes: 



ml 



[E - Nhoj -(P-(N- 2r)hk) 2 /2M]a r = — ^(y/r{N - r + l)a r 

8s 



+ Na r + y/(r + l){N-r)a r+1 ), 



(24) 



[E - Nhuj - (p - 2lhkf /2M]ai 



— — — a;_i + 2ai + ai +1 ). 
loe 



(25) 



This approximation preserves the nature of conjunctions between the neighboring ground level amplitudes: l — > l ± 1. 
Therefore, we would expect that this approximation does not cause qualitative changes for ground level stationary 
states. But it, as a matter of fact, strongly depresses the excited-levels, created from the excited-level atòmic states. 
The probability amplitudes, calculated on the basis of Eqs. (|2ÍJ) and (|2q), are plotted in Figures 5a and 5b. All the 
parameters are like in Figures 4a and 4b and have to fulfill the condition ( ^3| ) (|e| /(y/Nílo/A) — 58): Note that the 
chosen value |e| is even larger than the natural linewidth of envisioned optical transition. 

Unfortunately, even approximated Eq. ( |2~i| ) (or (|2^) ) can not be solved analytically due to (2rhk) 2 /2M terms, 
quadric relative to variable r (or ï). To proceed, we make a new supplementary approximation, regarded to as 
a Tavis-Cummings-type approximation |Ï3|| . Presuming the a r (;)-amplitude as a slowly varying function of r (l) , 
a permissible step for N 3> 1 and if other parameters are out of Bragg and Doppleron resonance conditions, this 
approximation envisages the variable r (/) as a continuous one. It should be expected, of course, to be a good 
approximation for a r m , when the momentum distribution is rather dissipated. 

Expanding the a r ±\ in a Taylor series up to the second order, 



G 



&r±l — fl r i 



da r 1 d 2 a r 
dr 2 dr 2 ' 



(26) 



and substituting this expression into the Eq. (|2J) (CPW-case) we get a second-order ordinary differential equation 
for the seeking probability amplitudes and energy vàlues 



d 2 A(x) 
dx 2 



4 



x 2hk r 



dA x 
dx 



4E r 



E — Nhuj 



A(x) = 0, 



(27) 



where x = (P/2hk) + r - (N/2), U 



-(MÏq/Ae) and the new amplitude A(x) is introduced by means of a relation 



A(x) = y/r\(N -r)\ 



The respective equation in SW-case is 

d 2 a(x) 
— + 

dx 2 



4E r 



E-Nhw 



a(x) = 0, 



(28) 



(29) 



where x = (P/2hk) + l, U being the same as in the previous case. The classical analog of U is the average potential 
energy of atòmic cm. in the light field. 

Note, that the resulting Eqs. (^) and (|2^) are no longer equivalent to the preliminary Eqs.(^) and ((|^) and some 
specific quantum features of the atòmic states may be lost (see in particular the last paragraph of Section IV). 

Let us first consider the more simple SW-case (Eq. (|29|)). We restrict ourselves with the so called red-detuning case 
when £ -< (U >- 0). Then the Eq^^]) is the equation for parabòlic cylinder function D v (y) with a real argument. 
Its two independent solutions are UA\ 



«i = D v (y), 



a,2 = D v {-y), 



where 



y = 2 



i/j 



x = 2 



1/4 



(30) 



(31) 



1 E - NHlü + U 

The requirement for D u (±y) to be limited for any vàlues of y (variable /) leads to a condition v = n 
Then in accordance with (p3), the energy gets only discrete vàlues and is determined by the formula 



E n = Nhw -U+ \fE r U (n+- 



(32) 
0,1,2,.... 

(33) 



Simultaneously both relevant parabòlic cylinder functions turn to a Hermite polynomial H n (y/y/2) and as a sequence 



y/2 n n\ 



H n V2 



1/4 



(2^ + ^ eXP 



u \2hk ) 



(34) 



where C n is a normalizing constant. It is worth to note that the size of quantization in (|3^) exactly copies the known 
frequency of time oscillations in the problem of atòmic diffraction in the classical standing-wave field, taking into 
account the cm. motion along the field direction. 

In CPW-case, in distinction with SW-one, we arrive at 



E n = Nhw + 



P 2 



1/N 2 



2M 1/N 2 + E r /U 



U 



h 2 ül 
16e 2 



UE r n 



(35) 



and 



exp 



i (ÍL - r ) 

JV V 2 ' ) 

y/rl(N -r)\ 



Hm 



y(r) ' 

V2, 





y 2 {r)~ 


j exp 


4 



(36) 
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where 



1 £ r \ 1/4 / P N 1 1 P 



v{r)=2 \YP + lT) {2M +r -^'l/N^ + E r /uJP2hk)- (3?) 



CP W energy spectrum (|35j) has two distinguishing features in respect to SW spectrum (|33|) . The first is, that the 
total momentum P enters into the expression of energy. This means that the atòmic c.m. motion is not bounded in 
stationary states in general (P ^ 0), but keeps some one-sided mean velocity. States become bounded only in the 
limit N 2 — ► oo, that is, in classical field limit. As a scale of this limit is claimed the ratio U/E r . Taking into account 
the definition of U (just after Eq.(p7j)), the mentioned condition may be written in the form 

ml , 

N>> im> w 

provided, of course, the condition (|23|). For V = 10~ 3 cto 3 cavity the condition ( j38[) is satisfied at N ^ 10 for sodium 
àtoms and at N )p 50 for cesium àtoms. 

The second difference of ([35]) relative to ( |33|) is the size of quantization, emerging by the additional term h 2 ÇL^/l&£ 2 
in square root expression. It has negligible influence in the viewed off-resonance approximation. 

Evident differences occur in probability amplitudes. It is convenient to image the expression in ( |36| ) as a product of 
two parts; first, including the Hermitc polynomial and second, the ratio of exponential to the square root of factorials. 
First part coincides with the (|||), but is shifted relative it by 

1 1 " 09) 



Í/N 2 + E r /U N 2 2hk 



This shift has the same coefficient at P/2hk as the kinetic energy term in (|35|), and is not small in quantum òptics 
range (N ~ 10 for off-resonance approximation). The second part, i.e. the additional ratio, is a symmetric function 
relative to point r = N/2, the form of which is presented in Fig.|| (for some vàlues of N). As is seen, it can be viewed 
as unity in central region and decreasing near the edges of the definition range. Hence, the amplitude distribution in 
CPW-case, being product of two symmetric by itself, but shifted with respect to each other functions, is not, already, 
a symmetric one, as it was in CPW-case. It is symmetric only if P = 0, but even in this case, the CPW distribution 
is somewhat narrower than the SW distribution. In analogy with energy spectrum, all differences disappear in the 
limit A 2 — > oo, where the counterpropagating light waves appear to a far detuned two-level atom as a sinusoidal 
potential, just similar to the standing light wave. 

Figures 7a and 7b show the probability distributions calculated by means of formulas ( |36| ) and ( |34|) respectively. All 
the parameters are like in Figures 5a and 5b. Comparison of Fig.7a with Fig.7b shows a small, almost imperceptible 
narrowing in each CPW-" mountain range" . The mutual comparison of Figures 7a and 7b with Figures 5a and 5b 
corroborates the Tavis-Cummings-type approximation to be an acceptable tool in the range of bounded states. But 
it totally puts out of consideration all free or quasi-free states. In Fig.8a we present a similar to Fig.7a distribution, 
calculated for conditions more precisely satisfying Tavis-Cummings-type approximation (A = 24). Its coincidence 
with the original, see Fig.8b, is of course better. However, it is to be memorized that the rising of N drifts the theory 
to the semiclassical domain of interactions. 



IV. FORMATION OF STATIONARY STATES 



Let us return to the theme of stationary states in the cases under consideration. Our objective is to clarify the ideas 
about the quantum contents of interaction mechanisms in SW and CPW-cases, to indicate the differences between 
them and already from this point, to understand the formation of stationary states from initially free states. To 
present them, we will first consider the SW-case. The variation of atòmic c.m. motion in this case can be imagined 
as a diffraction of atòmic matter waves on "frozen" immovable object: the standing wave grating. This is a fully 
classical-like phenomenon, provided by the matter-wave representation for atòmic c.m. motion. The picture, that the 
scattering is a result of absorption and reemission of SW-photons should be regarded as a defective, since it does not 
answer to the question how the atòmic momentum is changing due to these zero-momentum photons. This obstacle 
is overcame in the theory by the fact that the probability amplitude distribution ai with l = 0, ±1, ±2, ... , as well as 
the amplitude distribution bi, represent only the atòmic states and do not concern to photons. 

To realize the adiabatic evolution of the system from initially free state to the stationary coupling one, we assume 
the parameter Oq of the interaction to be a slow function of time, replace E by ih{d/dt) + E in ( p5| ) (in off-resonance 
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approximation) and view the az-amplitudes and the parameter E as a slowly varying in time quantities. After a 
phase transformation 



ai = ai exp 



Nhu -(p+ 2lhk) 2 /2M - 



8e 



we arrive at the equation 



Nül 

8e 



(cos(^_!/ 2 í) Im(aza:*_i) + cos(i/ í+1 / 2 í) Im(a;;a* +1 ) + 



(40) 



(41) 



sin(^_ 1( / 2 í) Re(aia*_ 1 ) + sin(t/, +1/2 í) Re(ona* l+1 )) 



for density matrix elements, which are more preferable for the putted objective. Here v\ = 2hk{p + 2lhk)/M is the 
analog of Doppler frequency shift. 

The equation (^ÏJ) is, of course, too complicated to be inspected in general form, but we may consider the initial 
stage of interaction, taking in addition ai(t — 0) = <5u - Then ol\±\ are fully imaginary with Imam >- (for e -< 0) 
and the relative rate of depletion of the state with / = Iq — 1 and l = Iq + 1 is given by relative vàlues of cos(f; Q _i/ 2 í) 
and cos(z/; 0+1 / 2 í) respectively. If p + 2lç,hk ^ 0, then cos(z%_i/ 2 í) ^ cos(z^ 0+1 / 2 í) and less modulo momentum 
states are being filled more rapidly than more modulo momentum states. As a result, the mean atòmic momentum 
decreases gradually, gathcring for asymptotically long times about the zero value. In attained stationary state, the 
mean atòmic momentum approaches zero. Note, that this is a stimulated analog of well known Doppler cooling 
(slowing) of red-detuned àtoms in the standing wave. The case of CPW we will start (and mainly restrict ourselves) 
by the immovable atom model (about commutative set of operators see in Appendix), since it g ives an op port unity 
to present the additional mechanism for stationary state formation. Putting M = oo in Eqs. (15a) and (15b) and 
substituting the second into the first we readily arrive at equation 



d . ,2 



ml 



A[E + he- NHuj] 
for density matrix elements, where 

a r = a r exp 



^Jr{N -r +l)lr[\{a r a r _ l ) + yj(r + Í)(N - r) lm(a r a* r+1 ) 



-- [ E-Nhw- 

h 



h 2 ülN 



8[E + he - Nhu] 



(42) 



(43) 



Here, in contrast to SW-case, each atòmic cm. state is single-bonded with the pair of photon numbers in CPW. The 
mathematical response of this connection is that the a r (and b r ) amplitudes are the representation of both atòmic 
cm. and photonic states simultaneously. The difference between the rates of r — ► r — 1 and r — > r + 1 transitions 
now is determined by square roots y r (N — r + 1) and y/ (r + 1)(N — r). The former term represents the process of 
reemission of a photon from the first travelling wave, containing preliminary r photons, into the second one, with 
respectively N — r photons, concomitantly altering the atòmic momentum in 2hk. The second term represents the 
inverse process, i.e. the reemission of a photon from the second wave into the first wave, accompanying it with a 
2hk atòmic momentum change. The physical reason of differences between the mentioned rates is solely the bosonic 
nature of photons. Direct comparison shows the dominance of reemission of photons from "more photon" wave into 
the "less photon" wave, which finally results in equalization of mean photon numbers in both counterpropagating 
travelling waves ]Ï5[ ], Hence, in attained stationary states the mean field momentum approaches zero. 

Thus, both analyzed cases possess diverse mechanisms of approaching to stationary state, but the final results are 
similar in the sense, that since one of the coupling subsystems, atom or field, is "immovable" or "frozen", the mean 
momentum of the other subsystem unavoidably approaches zero. Speaking in fancy form, the "immovable" part of 
the system takes on itself the "surplus" of the momentum from the other subsystem. A clear and expected result. 
In CPW-case with a moving atom, however, this result is not the general case. The reason is the conservation 
of system's momentum. In this case the system evolves under action of both mechanisms. First of them, the 
diffraction, tends to suppress the mean atòmic momentum, whilst second one, the bosonic nature, in contrary, tends 
to suppress the difference between the mean numbers of counterpropagating travelling waves. Since they are not 
possible simultaneously (except the case P = 0), in the course of time they mature the stationary states with definite, 
asymmetric in general, statistical distributions, where these tendències are mutually compatible. The role of bosonic 
mechanism however is essential in the quantum-optic domain (N =<! 10) and gradually diminishes out of it. 

To conclude this section, we would like to make a notice about one unavoidable peculiarity of Tavis-Cummings-type 
approximation, with respect to adiabaticity of interaction. Taking the photon momentum as a continuous variable, 
we have suppressed the size of its quantum up to zero (as small as possible) . Thereby we make feasible the arbitrary 
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small variations in the energy of atòmic cm. motion. As a consequence, any as small as possible variation of the field 
intensity (more concrete of parameter ÍIq) should cause changes in atòmic cm. motion and breaks dow n the a diab atic 
appr ox imat ion. Because of this, the approximate solutions (|34|) and (|36|), in contrast to exact Eqs. (15a), (15b) or 
( 18a ), ( 18b ), cannot be adiabatically related to the free system's state. Continuous quantum transition in atòmic 
cm. motion completely "washes out" the information about the initial state and its energy. Really, the switching 
off of interaction (ílo — * 0) in Tavis Cummings-type solutions leads to only one -E(ílo) = Nhoj value of energy and 
to one corresponding probability amplitudc with r = (N/2) + (P/2hk) photons in the first travelling wave and the 
rest N — r = (N/2) — (P/2hk) photons in the second travelling wave. It is necessary, however, to clearly realize that 
the parameters n\, n 2 and p preserve their physical essence as parameters of initial state; they are determined by the 
relations (0), ( |ï(ï| ) and (^3|) (in CPW-case) and are í^o-independent. Simply, the Tavis-Cummings-type approximation 
interrupts any, even adiabatic, connection of the stationary state with the initial one. The parameters ni,n 2 and p 
can be used in this approximation as externally (additionally) given ones, and in case of necessity can be corroborated 
by means of exact equations. 



V. BRAGG AND DOPPLERON RESONANCES 

These resonances, Bragg and Doppleron, have been identified for an atom moving in a standing wave | fÏ6f and have 
got implementations in atom òptics and làser cooling |Ï7| . Here we would consider only the limit of small intensities, 
focusing on the role of multiphoton recoil energy. 

Somc prcliminary words, for reminding, how these resonances are being extracted from equations. Thinking of 
the limit fio —> as a lcading one to the initial state, one should put a corresponding condition on the initial state 
amplitudes. For example, if ri\ and n 2 photons in CPW and the atom is on the ground level, then a r — <5 r . ni , b r = 0. 



As a consequence, we get from (15a) the value of the free system energy, 

E(Q = 0) = p 2 /2M + (m + n 2 )hw. (44) 

Here the atòmic momentum p is arbitrary. When in contrary, the limit ÍIq — s- is thought of as a transition to the 
final free state, we may only confirm the equality of energy to its (^J) free value. It follows from the conservation of 
energy. Also taking into account the relations (0) and (0), from Eqs. ( |l5a| ), ( |l5b| ) we directly arrive at 

(ni - r) (p + (-m + r)hk)a r = 0, (45) 



(lü-üj - (2n x - 2r - 1) kp/M - (2n x - 2r - l) 2 hk 2 /2M)b r = 0. (46) 

The solution with r — n\ and respectively a ril = 1, a n= ^ ní = 0, b r = 0, returns the system into the initial state. The 
next solution with 

p = (-ni + r) hk (47) 
gives a sole nonzero value for a r=ni + p /rife ground-level amplitude. Corresponding to this case final atòmic momentum 



(see the wave function (11)) pj = P + (N — 2r) hk = p + 2 (n\ — r) hk — — p, which is the condition for Bragg 
scattering. 

The third possible solution of (|ï|) and (|46| ) oceurs at 

lü - lüq = (2ri! - 2r - 1) kp/M + (2 ni -2r- l f hk 2 /2M (48) 

and has a sole nonzero solution for the corresponding excited-level amplitude b r . This is the case of Doppleron- 
resonance scattering, when the atom leaves the zone of interaction on excited energy level p8| . 

Let us stop our attention on Doppleron-resonance condition (|4Sh in order to examine the behavior of final momentum 
p f versus the initial momentum p. To this end the condition (Bq) may be viewed as a dependance of r on p and then 
put into the excited-atom final momentum expressionp/ = P+(N — 1 — 2r) hk = p+(2n\ — 1 — 2r) hk. Computation 
gives 



PI = -pv / l + 2Mh(Lü-Lü ) /p 2 , (49) 

where the atòmic initial momentum satisfies the Doppleron resonance condition ( |4"§| ) . For a fixed frequeney value uj the 
allowed vàlues of initial momentum p constitute a discrete set, corresponding to the integer vàlues of (2ni — 2r ~ 1) 
in @. 
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Graphs of (E9) are plotted in Fig]|. For comparison we also plot the graphs without the recoil energy term (the last 
one in (^)). As we can see, the presence of recoil energy dramatically changes the character of Pf(p) dependance. 
The first distinctive feature is in splitting of possible vàlues of pj for every p into two branches, herewith converging 
the high-order resonances (p 2 /2M ^ h \u> — ujq\) to the Bragg resonances (pf = —p) in the additional branch. The 
second difference is the lack of resonances in small atòmic momentum range (p 2 /2M -< h\u> — u>o\) and connected 
with it, the convergence of p/ to zero at the limiting point with p 2 /2M = h\u> — ujq\ . It means that the inherent 
possibilities of the Doppleron-cooling really are more than it follows from the theory, not including the recoil energy 
term. 



VI. CONCLUSIONS 



We analyze the stationary states of JCM with SW and CPW fields, treating the external, as well as the infernal, 
atòmic coordinates quantum-mechanically. The inspection shows essential deviations between SW- and CPW-cases 
in quantum òptics range (N -< 20), which first of all have qualitative character. As a such matter can be noted 
the deterioration of the concept of quasimomentum, the asymmetry and narrowing of momentum distribution, the 
coincidence of photonic and atòmic momentum distributions in CPW-case. 

Since the physical reason of the mentioned differences is the existence of the momentum conservation, it can be 
concluded that the modern cavity QED an the theory of measurement, being developed on the base of quantizcd 
SW-cavities, in fact have only dealt with the conservation of energy. The cavity QED and related problems for ring 
cavities, where the momentum conservation exists too and has to alternate the state-pattern in " atom+field" system, 
need a performance, yet. Moreover, the experimental confirmation just of CPW-case regularities, stationary or not, 
may be adopted as the direct nonmediated evidence for the photon momentum to be quantized. 

Reexamination of kinematic resonances shows that due to recoil energy term the Doppleron-resonance splits into 
two branches and , which is more important, can not take place in the range of small incident kinetic energies (where 
H -< ^/2Mh\Lü~uj \). 

This work was supported by ISTC Grant A - 215 - 99. 



VII. APPENDIX 

The case of immovable atom (p 2 /2M = 0) is a separate one for stationary state problem and its analysis is not 

, ( |l5b| ). The reason for such a separation is that the P and T operators loose the 
inherent physical contents and must be omitted. Instead of them enters a new, commutative with H and N photonic 
operator 

9 = {a^ai + a^a 2 — a^a 2 — aJai)/2iV, (Al) 

which has the eigenstates 6 m = 1 — m/N, m = 0, 1, 2, • ■ •, N. 

It is referred to as an operator of interference, in the sense that its eigenvalues play the same role in the quantized- 
field picture as the term cos 2 kz in the classical field picture of interaction |Ï6[. 
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FIG. 1. Dispersion curves of the atòmic cm. quantum motion in countcrpropagating waves (a) and standing wave (b) 
quantized cavity fields with iV = 6 photons. The chosen parameters are he/E r = —250, Ç = hQ.o/2 3ll2 E r = 1.75. The latter 
corresponds to 3Si/2 — 3P 3 / 2 optical transition of sodium atom located in a V = 10~ 3 cm 3 quantum cavity. n = 1, 2, ... 
on curves labels the number of energy zone. The lowest states, up to n = 4, are shown. CPW-curves (a), in distinction to 
SW-ones, aren't periòdic and, as a result, the concept of quasimomentum can not be introduced for them. The energy distance 
between the n = 1 and n = 2 curves, that is the width of first energy band, is much wider in CPW-case than in SW-case. 

FIG. 2. Dispersion curves of the atòmic cm. quantum motion in counterpropagating waves (a) and standing wave (b) 
quantized cavity fields with N = 12 photons. The other parameters are like in Fig.l. The wider widths of energy bands in 
CPW-case still is aparent. The second band between n — 2 and n — 3 zones is clearly seen in CPW-case, whilst in SW-case 
it is not. 

FIG. 3. Distribution of the atòmic cm. momentum states over the energy spectrum. The probabilities are presented on 
vertical axes. The momenta are on the frontal axes and are presented in p r = 2hk units. The energies are presented without 
photonic field energy in E r = (2fifc) 2 /2M recoil energy units. The total momentum P (quasimomentum p) of "atom +cavity 
field" system is chosen zero. The other parameters are like in Fig.l. The vàlues of energy are quantized and depicted by dark 
lines on both graphs. The spectrums of momenta are quantized too, but for sake of simplicity of exhibition we connect the 
neighboring points in both momentum and energy directions. The lower-laying energy levels are created from ground-level 
free atòmic states, and the upper-laying ones are from excited-level free atòmic states. The energy spectrum of CPW-case (a) 
is not truncated, whilst the spectrum of SW-case (b) in reality is infinite and we exhibit only some truncated parts for both, 
lower-laying and upper-laying famílies. Finally the distributions in a form of expanding hills pertain to quasi-bounded states; 
the others are localized near the zero-momentum vàlues and represent the quasi-free states. 
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FIG. 4. Distribution of the atòmic cm. momcntum states over the energy spectrum for N — 12 photons in cavities. The 
parameters are like in the previous figure. The momentum distributions in bounded states are somewhat wider than in previous 
N = 6 case. The distinction of CPW-case (a) from SW-case (b) is essential yet. 

FIG. 5. The same momentum distribution as in Fig.4, but calculated in off-resonance approximation. Excited-state energy 
levels are escaped, of course, while the ground-state energy levels and c.m. momentum distributions practically coincide with 
the corresponding results in Fig.4. 

FIG. 6. The function-factor, which determines the form distinction of c.m. momentum distribution in CPW-case with respect 
to SW-case (in Tavis-Cummings approximation) , as function of number r of photons in one of the counterpropagating waves. 

FIG. 7. Atòmic c.m. momentum distribution in Tavis-Cummings approximation for CPW (a) - and SW (b) -cases. The 
parameters are taken from Fig.5. The continuous spectrum absent in this approximation and the difference between CPW- 
and SW-cases is notable only for upper-levels. 

FIG. 8. Momentum distribution for atòmic c.m. quantum motion in off-resonance and Tavis-Cummings approximations with 
N — 24 photons in cavities. As was expected the coincidence is better (with the lower bounded states) than in the previous 
N = 12 case, but is not appropriate for quantitative calculations. 

FIG. 9. Final atòmic momentum as a function of initial (incident) momentum for Doppleron-resonance conditions with 
(sòlid line) and without (dashed line) recoil energy term. The graphs are depicted as continuous ones (instead of discrete point 
sequences) only for the sake of simplicity. 
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